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ABSTRACT 



We investigate the physics of elementary excitations for the so called relativistic scalar 
plasma system. Following the standard many-body procedure we have obtained the RPA 
^. equations for this model by linearizing the TDHFB equations of motion around equilibrium 

r-| ! and shown that these oscillation modes give one-meson and two-fermion state of the theory. 

The resulting equations have a closed solution, from which one can examinate the spectrum 
K^ ; of excitation modes. In partticular, our results indicate existence of bound state for certain 

\-{ ' region of phase diagram. 
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I. INTRODUCTION 

In a previous work [^ (hereafter referred to as I) we have presented a framework to in- 
vestigate the initial-value problem in the context of interacting fermion-scalar field theories. 
The method allows one to obtain a set of self-consistent equations for the expectation values 
of linear and bilinear forms of field operators. These dynamical equations acquire the struc- 
ture of kinetic type, where the lowest order of the approximation corresponds to the usual 
gaussian mean-field approximation (coUisionless). As application, we have implemented a 
zero-order calculation within the simplest context of relativistic scalar plasma system. We 
have shown that the usual prescription of renormalization can also apply to these nonop- 
erturbative calculation. In particular, we have obtained a finite expression for the energy 
density and the numerical results suggested that the system presents always a single stable 
minimum . 

In continuation of I we will report in this paper a particular application of the renormal- 
ized nonlinear obtained the previous publication. We follow here a recent work by Kerman 
and Lin [0, |^ in order to investigate the near equilibrium dynamics around the stationary 
solution. We shall show that one-meson and two- ( quasi )fermion physics can be studied from 
the linear approximation of the mean-field equations. In particular, one can solve these 
equations in a closed form and find scattering amplitude as well as the conditions for the 
two-fermion bound state. 

For completeness and notational purpose, we repeate here a few key equations of I. A 
summary of derivation for these equations is shown in Appendix A. For the scalar plasma 
model, the dynamics are governed by the hamiltonian 



H = n 



L 



Ti = —ipli'j.d — irijip — gip(t)ip H 



|^n^ + (i + z)|a0p + /.V 



(1.1) 



[We use the notation: /^ = j d^x] where ■j/' is a spin-i field while is a scalar field. The 
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parameters m and /i are, respectively, the mass of fermion and scalar particles and g is the 
coupling constant. The last term of this expression is the counterterms necessary to remove 
the infinities occurring later [jl|, Q, 

4-Hc = ^0 + |^0^ + §0-V^0^ , (1.2) 

where the coeffcients A, Sfi^, C, D, and Z are infinity constants to be definided later. 

In order to study the dynamics of the system, we focus on the selected set of observables, 
which are related to the expectation values of linear and bilinear forms of field operators, 
referred to as gaussian variables. The time evolution of these quantities obeys the Heisenberg 
equation of motion, 

i{d)=Tr^40,H]J^ (1.3) 

where O could be 4>{x), (f){x)(f){x), ijj{x)ip{x) and etc., and F is the state of the system in the 
Heisenberg picture. As approximation we replace the full density JF by a truncated ansatz 
jFo(t). By construction, jFg reproduces the corresponding JF averages for linear or bilinear 
field operators [see Eqs.(41) of [|l|]. In particular, we have used a formulation appropriate 
for the many-body problem, so that JFq can be written in the momentum basis as & 0] 



B tF 



J-Q — J-qJ-q 



T 



1 / V. ^''^'^ 







rr^_^:^ (1.4) 

J-„-^ 1 + z/„ I 1 + z/„ i ^ ^ 



+ z/p V 1 + ^] 

,(A)^(A)t (A) , ., „Wx ,(A)^(A)t, 



•^0^ = nKi<i<i + (i--P<HT] , (1-5) 

k,s,A 

where a (a^^) and P (P^^) stand for Bogoliubov quasi-particle annhihilation (creation) op- 
erators for fermion and boson respectively. In doing so, the gaussian variables are now 
represented by the Bogoliubov parameters [see ( |A.3D -( |A3D ] and its equations of motion, can 



be obtained directly using ( |1.1|) - (|1.5|) . The resulting expressions are: 



where 



0k = fi'(0)^— sin7k (1.6) 

Kq 

7k = 2ko-2gm{(f))—[l cot2v9kCOS7k (1.7) 

ko \ m J 

W = jj|^(n> (1,8) 



£w^-^W'-(|^ !/.(-) + /.(-)! (1.9) 



^^"^ = (2^/'''P2(?T?F5 <1'1°) 



J3, " 



^i('") = /'"^ (k^ + m2)i/2 '=°"^^'- (!■") 

|k| 
/2(m) = /rf^'k-^^^^— ^sin2v9kCOS7k . (1.12) 



The mean-field energy density, on the other hand, reads as 

{H) 1 



V V 



TrHTo 



-2 / ko cos2^k + 2^(0) [him) + him)] 



The equations above describe the real-time evolution for the relativistic scalar plasma 
system in the gaussian mean-field approximation. They are nonlinear time-dependent field 
equations. Therefore, a closed solution is not easily constructed. Here, we will consider the 
equations in the equilibrium situation and the small oscillation regime. In these cases, a 
closed solution can be found allowing us to examine diverse properties of the theory. The 
structure of this paper is as follows. In Sec. II we shall derive the RPA equation for this 
model by considering near equilibrium dynamics about the equilibrium solutions obtained 
in I. Section III shows an analytical solutions for these RPA equations using a well know 
procedure of the scattering theory. In Section IV we discuss the question of renormalization 
within the context of scattering amplitudes and examinate the possibility of existence of 
bound state soultions. 

II. NEAR EQUILIBRIUM DYNAMICS 

As a function of the Bogoliubov parameters, the energy density is a function of these 
variables. A minimum of (|1.13|) corresponds to the ground state of the system and small 



amplitude motion about the minima is described by the linearized approximation of the 
gaussian equations motion ( |1.6| )-( ]L9| ), yielding a set of harmonic oscillators ^. The normal 
modes and the eigenvectors are the RPA solutions. In the field theoretical context these can 
be seen as the energy and the wavefunctions of quantum particles. In this section we derive 
the RPA equations of the model and next section will discussion solutions for this problem. 
Let us consider first the static problem (see section Vl-b of I). Recalling (|1.6|)-(pr9|), the 
stationary conditions require 7k = 0k = /^p = % = {<P) = (H) = 0. In Section IV-b of I 
we have discussed in detail the renormalization conditions and the solutions and for this 
set of equations. In particular, the following choice for the coefficients of He satifies the 
self-consistency condition [|l], || 



D 

C 
A 



±A87Tg^L{m) , 

T247i^g^L{m)G{fi) T 167r^^G'(0) ± 247im^ g^ L{m) 

^ASnmg^ L{m) , 

±24:7img^L{m)G{fi) ± IGnmgG^m) , 



with 



L(m) = 



1 



/k 2k2(k2 + m2)V2 • 
In doing so, the finite static equations of the system can be regrouped as follows: 



(2.1) 
(2.2) 
(2.3) 
(2.4) 



(2.5) 



sin7k|eq = 



cot 2ipk 



(n)U = o 



■^fJ'^{<P)\c<i-9m'^ 



(k + mm) 



Iklfm — m) 



ln(^ 



(2.6) 
(2.7) 
(2.8) 
(2.9) 



where m is the effective fermion mass, 

m = m-g{<f))l^ . (2.10) 

The equation (2^)-(2^) can be solved numerically for any given values of /i and g in unit 



of m. 

Investigation of the near equilibrium motion proceeds by considering flutuactions around 
the stationary point 
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<y^k = ¥^1°^ + ^<y^k 



7k = 7^°^ + '^7k 

(2.11) 
(0) = (0)(°) + 5(0) 

(n) = (n)W + 5(n) , 

where ^'^\l^^\ (0)^°^ and (n)W satisfy Q-Q and the quantities 5(^k, 57k, 5{(p)l^ and 
(5(n)|^q will be assumed small in our approximation. Next step is to expand ( |1.6| )-(|T^) up 
to first order in these fiuctuations to yield 



50k = ^(0)(°)^57k (2.12) 



^(0)(°)|k|57k = -^ko{k'' + My^^-2g\k\koS{^) (2.13) 



m = (T^^W (2.14) 

Eliminating the momenta 57k and 5(11) we might rewrite ( p.l2|) -( p.l5|) in a more compact 
second order equations: 



5(^k = -4a;'5v?k-2^|k|5(0) (2.16) 



{l + Z)5{4)) = - (/x^^ + S) (5(0) + IGirg J h{k)5ip^ (2.17) 

where we use notations 

h{k) = B (2.18) 



with /cq = \^^^+m? and 

S ^ <5/x2 + ^4vrG(/i) + C(0)(°) + ^(0)(°)' . (2.19) 

Thus, the small oscillation dynamics is described by coupled linear oscillator equations as 
usual in RPA treatment. In particular, these modes decouple when g = yielding two simple 
oscillator equations. 

A solution to this problem involves determining the normal modes of small oscillation 
and their frequencies. On the other hand, earlier studies have been demonstrated that these 
elementary excitations can be interpreted as quantum particles. In our case here, (5(/9k can 
be seen as two-fermion spinless wavefunction while 6{(f)) gives an one-meson physics of the 
system in this regime. Furthermore, k stands for the relative momuntum of two fermions 
and its total momentum is constant of motion for this context of uniform system, therefore 
there is no explicit dynamics involved. Notice also that for scalar setor, the particles do not 
depend on the momentum in this context. 

III. RPA EQUATIONS AS A SCATTERING PROBLEM 

In the last section, we obtained the linear approximation for the gaussian equations of 
motion. They describe the physics of elementary excitations of the system. We shall show in 
this and next section that this coupled linear oscillator equations can be solved analytically 
giving the wavefunctions and the spectrum of excitations. 

Let us first Fourier transform the wavefunctions to the energy representation. 



(3.1) 
S{(/)){t) = f duj5{<P){uj)e'^' 



where Sif]^ and (5(0) are now energy- dependent amplitudes. Substituting now (|3.1|) into 
(|1|)-(|13) we have 



u'^ - Au"^) 5^],{u) = 2g\k\6{(f)){uj) 



(3.2) 



(cu^-/i^ + Zo;2-S)5((/))(cj) 



-I6ng / /i(k)(5v?k(w) 

Jk 



(3.3) 



Since the oscillation amplitudes play the roles of wavefunctions of quantum particles, it is 
then more convenient to treat this system as a coupled channel scattering problem with 
appropriate boundary conditions 0. 

Next, we will discuss the scattering process where the source is a two-fermion wave. In 
this case, we can solve fB.3|) as follows: 



^(0)H 



— IQirg 



h(k')6ip]ii{uj) 



(3.4) 



u;^ — /i^ + Zlu'^ — T. + irj 7k' 
In this expression we have included the boundary condition of there is no incident wave of 



(5(0). Substituting this into ( |3.2| ) and rewriting it in terms of a new variable 

^k = ykoSfk 
one finds 



(^^ - ^kl) ^, 



-327Tg^ 



-.^v 



(3.5) 



(3.6) 



u^ - n^ + Zuj^ ~ j: + 17] ^^ A' y^ 

where the Green's Function includes the effects of coupling of ^P to (5(0). Notice also that 
the potential is separable in the sense that H] 



{k\V\k') = v{k)v{k') 



Ikl Ik'l 



ko Y Kq 
It is now convenient to rewrite it as an integral equations, 

*(k, q;u;) = a5{q-k) 



(3.7) 



+ 



1 



-327r^2 



[a;2 - 4ki + ie] [a;^ - ^-^ + Zuj^ - E + i7]] ' ' 7k 
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t;(k) / t;(k')vl/(k',q;a;), (3.8) 



where q is the relative momentum for two incident quasi-fermions and a is an overall phase 
factor. We have used the boundary condition of the outgoing wave condition (+^e) as solution 
of Eq. (|3.7| ), but we could have chosen e.g. the incoming wave condition (— ie) or Van Kampen 
wave condition or other conditions. 

The integral equation (|3.8| ) can be solved as usual. We first multiply the expression In 
order to solve we integrate this expression with respect to k, 

/ ^(k)^(k, q; u) = :^ (3.9) 



where 






^ (^) = / /YT^-2 . TJT-T ^3.10) 



Substituting this back result into (P78|) yields 



*(k,q;.,^„.(q-k, + -j-2_|!|I (3,1) 

with 

A+H = -^ {u' - ^2 + Zcu' - S) - J+H (3.12) 

We have, thus, found an analytical solution for the elastic channel of two-fermion scattering 
problem given by (|3.2|) - (|373|) . 

The special form of interacting potential in this case allows one also to get easily a closed 
expression for the scattering matrix 

Tik,k';u;) = v{k)^^vik') (3.13) 



with A+(u;) given by (p. 12 ). In summary, this section has discussed the solutions for RPA 



equations. These elementary excitations describe a coupled channel scattering problem and 
we have studies the particular case of two-fermion elastic process. For this simple interacting 
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potential, it is easy to obtain closed expression for the two-fermion wavefunction and the 
scattering matrix, where several dynamical behavior can be read off from A'^{uj). The 
remained problem is the divergent integral I^{uj) in ( |3.12|) which must be removed with the 
help of countertems. We shall show in the next section that besides the counterterms used 
static discussion also apply we will need a convenient choice of Z. 

IV. RENORMALIZATION AND BOUND STATE SOLUTION 



This section will investigate the conditions for the existence of bound states of Dirac 
spin-1/2 particles in a system of scalar plasma, utilizing the framework developed in the 
previous sections |T^. In this context, the standard procedure is through an analysis of the 
positions of poles of the scattering matrix. The equation ( 3.12| ), however, has a divergent 
integral. We shall show next that the divergent terms can be kept directly under control 
with the help of (2J-)-(2]l) and a convenient choice of Z, yielding a finite expression for 

Let us thus substitute the counterterms (P^.1D-(P^ into (|3.12|) , after some algebra one 

gets 



A+M 



ing' 



;i + A7Tg'L{m))uj' - /i^ + 167r^^G(0) - 247Tg^M'L{m 



r(uj) 



(4.1: 



with /^(cu) given by ( |3.1CI|) . Notice that in the interval of < u; < 2M the integral /^ is well 
defined and we can take e = 0. For u > 2m, on the other hand, the spectrum is continuum. 
The calculation is straightforward and we find: 



nu^) = Q'^F{u)-9iuj-Am')^ 



u? - Am^ 



where 



Q 



An 



A^ + — - 3m^ log — 
V 2 m 



and the finite terms are 



n^) 



uo^ - Qrfi^] log + ^ — tan-^ 



uo^ 



\2mJ 



UJ 



Am? — ll? 



(4.2) 



(4.3) 



2 ^ /1„^2 



< cj^ < 4m 
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(4.4) 



F{uj) = (m2_6o;2)log — +^ ^ log , ^ , ^ u? > An? 

^ ' \2mJ uj u-^u^-Am? 

(4.5) 



Comparing ( [4.1| ) and ( [4.2| ) one notes immediately that there is still a logarithmic divergence 



left. This will be canceled with the following choice of wavefunction renormalization [Q 

Z = Ang'^Lim) . (4.6) 

The resulting finite expression is 



,2 



A+(cu) = -4^ + F{^) - ^(^ - 4m2)— L^ - 4m'l . (4.7) 

Discussion of the problem consist now in solving the equation 

A(c<;) = 0. (4.8) 

Depending on the value of uj the system has different dynamical behavior. When u'^ < the 
system is unstable and for u;^ > we are in the scattering regime. Solution of interest here 
is in the interval of < tu^ < 4m^. In this case, the system may present a stabe bound state 
if one finds u^ such that A^{ujg). The fig. 1 illustrate A'^{uj) for several combinations of fi 
and g in unit of m. Notice that for the same value of g{= 1), A+(co') has a single(none) zero 
when fi/m < l.TdAlfi/m > 1.794). A natural way to interprete this result is that the meson 
mass determines the range of the Yukawa potential for a same coupling. When fi of the 
system is large, it is more dificult to the fermions to interact and consequently decreases the 
probability of forming a bound state, and vice-versa. One can, however, compensate with 
inreases of values oi g. In order to see this, we have examinate ( |4.7| )-( |0| ) for the parameter 
space fi/m versus g, the result is shown in the figure 2. 

In summary, this work has investigated the physics of elementary excitations for the so 
called relativistic scalar plasma model. To reach this goal we have derived RPA equations 
for this sistem by linearizing the TDHFB equations obtained in a previous work. In this 
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context, the amplitudes of excitations are identified with quantum particles and the resulting 
equations allow one to study scattering processes nonperturbatively. We have solved this 
RPA equation analytically using well known procedures of scattering theory where the scat- 
tering amplitude obtained has a simple form. We have also shown that the usual definitions 
of countertems can be applied to this resulting expression, from which relevant physics about 
the excitations of the system can be obtained. In particular, our results indicate existence 
of bound state for certain region of phase diagram. 

Apendix A: Mean-field kinetic equations 

In this apprendix we will review briefly the results obtained of I. The method used 
there was developed earlier in the context the nonrelativistic nuclear many-body dynamics 



by Nemes and de Toledo Piza |Tl|. More recently, our group has applied this technique 



to 0^ field theory |Q as well as to the Chiral Gross-Neveu model [Q. Our approximation 
focus on the time evolution of a selected set of observables, which in the case of scalar 
plasma system are theexpectation values of linear, 0(x), and bilinear field operators such as 
0(x)0(x) , il!{x)ip{x) , ip{x)ip{x) and etc. Because of further convenience, we work instead 
with expressions which are bilinear in the creation and annihilation parts of the fields. 
In the Heisenberg picture, 4>{x) are scalar spin-0 fields 



>^) = E 



{2Vpoy/' 



b^{t)e'P-^ + bl{t)e 



— ip.x 



(A.i: 



where &p(t) and bp(t) are boson creation and annihilation operators. For the spin-1/2 fields 
we have 









,(1)/+N„ik.x 



u,{\^,s)a'^^Ut)e'^-- + U2{\^,s)a'^nt)e 



(2)1 



-ik.x 



(A.2) 



,(1)1 



(2)/,\^ik.x 



^i(k, s)a\^:{t)e-^^-- + W2(k, s)a'Cl{t)e 



where a\^ [. (t) and aj^ {.{t) [a\^ {. [t] and ctk ^(t)] are fermion creation and annihilation operators 
associated with positive [negative]-energy solutions Mi(k, s) [M2(k, s)] of Dirac's equation. 
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In order to deal with condensate and pairing dynamics of the scalar field we define the 



Bogoliubov transformation as follows [[12 
dp{t) 



dlJt) 



cosh/tp + i^ 



— smh.Kp — i^ coshKp — t-j- 



sinliKp + i^ 



Vp 






(A.3) 



where dp is the shift boson operator defined by 



dp{t) = bp{t) - B{t)6pfl with Sp(t) = (6p(t)) = Trep [6p(t)^] 



(A.4) 



In the case of fermions we restrict ourselves for simplicity to Nambu transformation, which 
can be parametrized in a form that incorporates the unitarity constraints as 



" aS - 




«S 













cosyjk 



cos (y9k e *'"' sin (^^ 



—e^^^ sin y^k cos (p^ 



g*7k gi]2 y9j^ 







r ^(1) 1 


e ^'^^ sin (y9k 




«k,s 







«s 







(1) t 


cosv^k 







(A.5) 



The next step is to obtain the mean-field time evolution for the mean values of the 
gaussian observables in the context of the initial-value problem. In other words, we want 
the gaussian mean- field equations of motion for the parameters v^k(^) , 7k (^), Vpit): i^p{t)y 
for the quasi-particle occupation numbers i']i,\{t), z^p(t) and for (0) and (11). In Ref. |jl[| we 
obtained for the bosonic variables 



z>p = Tr,MPp,H]J^ 



ikp - e '^p {f]p + kpT]p) 



Tr[/5Mp,i7]^ 



l + 2z/o 



^(0(t)) = Tr^44>{t),H]J^ 



(A.6) 

(A.7) 
(A.8) 
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In the case of fermionic variables, we obtained 



i^KX = Trep([4f«g,^]-^) (A.9) 



[«V^k + 7k sm v^k cos v^kj e ^'' = ^^ ^^ ^^ . (A. 10) 



t'l 



k,s '^k,s 



Our implementation of the gaussian mean-field approximation consists in approximating 
JF by a truncated many-body density operator J^o{t) = J-'qJ-'q- The factorized form of the 
jFo(t) embodies what we refer to as the double mean field approximation. The subsystem 
densities JF^ and JFq are in fact unit trace gaussian densities, written in the form of an 
exponential of a bilinear, Hermitian expression in the creation and annihilation parts of 
the bosonic and of the fermionic fields, respectively. They can be, however, rewritten in 



diagonal form, (|1.4|) and (|1.5|) , when one uses Bogolyubov quasiboson and quasifermion 
operators. With these assumptions the equations ( |A.6| )-( |ATTo| ) of are closed now and they 
will determinate the time rateof change of gaussian variables. The results are shown in 

(PD-(PD. 
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Figure Captions 

Figure 1. The behavior of the function A{uj) as a function of energy in unit of m. 

Figure 2. Existence of bound state of two fermion as a function of parameters fi/m and g. 
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